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(57) Abstract 

Signals exchanged between a fixed 
infrastructure (1 10-1 16) and a mobile unit 
(160) give rise to M channel quality met- 
rics. The channel quality metrics are 
mapped, via predetermined relationships 
(501-503), into M corresponding time of 
arrival variances, which in turn are used 
to derive M-l time of arrival differential 
variances. A time of arrival differential 
weighting matrix including, in part, the 
time of arrival differential variances, is 
used to calculate a WLS solution, which 
solution is an estimate of a location of the 
mobile unit This procedure may be im- 
plemented using an infrastructure-based 
location processor (130) operating in con- 
junction with the mobile unit, or may be 
performed by either the location processor 
or mobile unit alone. 
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METHOD AND APPARATUS FOR COVARIANCE MATRIX ESTIMATION 
IN A WEIGHTED LEAST-SQUARES LOCATION SOLUTION 

10 Field of the Invention 

The present invention relates generally to location 
determination in wireless communication systems and, in 
particular, to the estimation of a covariance matrix for a 
15 weighted least-squares solution. 

Background of the Invention 

Wireless communication systems are known in which 
20 mobile units (e.g., in-car mobile or in-hand portable radios) 
wirelessly communicate with a fixed communication 
infrastructure comprising a plurality of geographically- 
diverse transceivers. In such systems, methods for 
determining location information for a given mobile unit are 
25 known. In particular, the well-known weighted least squares 
(WLS) solution can be used to determine location information 
as shown, for example, in U.S. Patent No. 5,416,712 issued to 
Geier et al. 

Succinctly stated, the WLS approach to location 
30 determination attempts to iteratively derive a location 
estimate for a mobile unit based, in part, on distance 
estimates between the mobile unit and fixed transmitters 
having known locations. Given that distance can be calculated 
as the product of velocity and time, the distance estimates 
35 (referred to as pseudo-ranges or PRs) are calculated in 
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practice by multiplying the propagation delays between the 
mobile unit and fixed transmitters with the speed of light 
Assuming ideally measured propagation delays, the location of 
the mobile unit can be calculated using the pseudo-ranges with 
5 little or no error. However, propagation delays are measured in 
practice using transmitted signals, which signals are subject 
to the effects of various error sources, e.g., noise, multipath 
interference, distortion, etc. The resultant errors in the delay 
measurements are translated into errors in the pseudo-ranges 

10 and, consequently, into error in the location estimate. 

In order to combat the presence of measurement errors, 
the WLS solution factors the reliability of the various 
measurements into the location estimation. That is, the WLS 
solution, when used to estimate location, places greater 

15 reliance on measurements having greater reliability, and 

discounts those measurements having less reliability. In this 
manner, the WLS approach offers a significant performance 
advantage over non-weighted techniques. In practice, however, 
the reliability of a signal cannot be directly measured, but 

20 must be described as the variance of a stochastic variable. 
Various methods for incorporating measurement 
reliability into a location determination context are known in 
the art. For example, in U.S. Patent No. 5,202,829 issued to 
Geier, Kalman filters are used to assess the "quality" of 

25 pseudo^ranges measured using GPS receivers located in a ship 
and accompanying tailbuoys. Additionally, in U.S. Patent No. 
5,436,632 issued to Sheynblat, a system is disclosed wherein 
redundant GPS receivers at known locations (reference 
stations and integrity monitors) are employed to provide 

30 corrections for pseudo-range measurements made by a mobile 
unit also equipped with a GPS receiver. Sheynblat discusses 
the use of a WLS solution in which "measurement error 
covariance" is determined based on receiver noise as 
represented by errors between a given reference station and 

35 its corresponding integrity monitor. 
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While Geier and Sheynblat incorporate reliability 
determinations to improve the accuracy of location estimates, 
such techniques are not readily adaptable to mobile and 
portable wireless communication environments. First, both 
5 Geier and Sheynblat require the use of GPS receivers. Such 
receivers add significant cost to mobile/portable equipment 
and prohibitively add to the size and complexity of such 
equipment, particularly portable radios. 

Additionally, as noted in '872 Geier, the use of Kalman 

10 filters requires significantly more computing power than a 
WLS solution. In portable units, computing power is often 
limited by size and battery life considerations, making the use 
of the computationally-expensive Kalman filter approach less 
attractive. Therefore, a need exists for a method which 

15 incorporates the advantages of a WLS location solution without 
the need for costly GPS equipment. In particular, what is 
needed is a technique for estimating a covariance matrix for 
use in a WLS solution. 

20 Brief Description of the Drawings 

FIG. 1 is a block diagram of a wireless communication 
system. 

25 FIG. 2 is a block diagram of a location processor. 

FIG. 3 is a flowchart illustrating operation of a mobile 
unit when the mobile unit cooperates with a location processor 
to determine a location for the mobile unit. 

30 

FIG. 4 is a flowchart illustrating operation of a location 
processor when the location processor cooperates with a 
mobile unit to determine a location of the mobile unit. 
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FIG. 5 is a graph illustrating examples of predetermined 
relationships used to determine time of arrival variances. 

FIG. 6 is a flowchart illustrating the determination of a 
5 location estimate using a weighted least squares solution. 

FIG. 7 is a flowchart illustrating operation of a fixed 
infrastructure when the fixed infrastructure determines a 
location of a mobile unit based on signals transmitted by the 
10 mobile unit. 

FIG. 8 is a block diagram of mobile unit. 

FIG. 9 is a flowchart illustrating operation of the mobile 
15 unit of FIG. 8 when the mobile unit independently determines 
its own location. 

Description of a Preferred Embodiment 

20 The present invention provides a method and apparatus 

for estimating a covariance matrix for use in weighted least 
squares location determinations. Generally, signals exchanged 
between a fixed infrastructure and a mobile unit give rise to M 
channel quality metrics (M > 3). The channel quality metrics 

25 are used to derive M corresponding time of arrival variances, 
which in turn are used to derive M-1 time of arrival 
differential variances. Predetermined relationships are used 
to map the channel quality metrics in to the time of arrival 
variances. A time of arrival differential weighting matrix-an 

30 estimate of a covariance matrix-comprising, in part, the time 
of arrival differential variances, is used to calculate a WLS 
solution, which solution is an estimate of a location of the 
mobile unit. This procedure may be implemented using an 
infrastructure-based location processor operating in 

35 conjunction with a mobile unit, or may be performed by either 
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the location processor or mobile unit alone. In this manner, a 
WLS solution can be advantageously employed in a location 
determination context without the need for GPS equipment. 

The present invention can be more readily described with 
5 reference to FIGS. 1-9. FIG. 1 is a block diagram of a wireless 
communication system 100, such as an IDEN™ system 
manufactured by Motorola, Inc. The wireless communication 
system 100 comprises geographically-diverse fixed 
transceivers 110-116 coupled to a location processor 130 via 
10 a network 140, and a plurality of mobile units 160 (only one 
shown). Although not shown in FIG. 1, each mobile unit 160 
includes memory and one or more digital computing devices 
used to store and execute software programs, as known in the 
art. 

15 The coverage area provided by each transceiver 110-116 

gives rise to a corresponding cell 120-126, represented by 
hexagons as shown. In practice, each transceiver 110-116 may 
actually comprise a bank of transceivers, although it is 
assumed hereinafter that each cell includes only a single 

20 transceiver for simplicity. While a mobile unit roams 

throughout the system 100, wireless communication services 
for the mobile unit are managed and provided by a serving site. 
Neighbor sites are those sites surrounding the current serving 
site. In the example shown in FIG. 1, the site identified by 

25 reference numeral 110 functions as the serving site for the 

mobile unit 160, and the sites identified by reference numerals 
111-116 correspond to neighbor sites. 

The transceivers 110-1 16-which in the preferred 
embodiment provide wireless channels according to a time- 

30 division multiplexed (TDM) format-and the location processor 
130 collectively are referred to as a fixed infrastructure. 
Being fixed, the precise location of each transceiver 110-116 
is known with high accuracy as determined by surveying or 
other suitable measurement techniques. 
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As described below, the present invention relies on 
channel quality metrics derived from signals received or 
transmitted by a mobile unit. In the preferred embodiment, 
each of these signals comprise known data symbols unique to 
5 each transceiver. For example, in an IDEN™ system, each 

transceiver will transmit its known data symbols whenever an 
unused time slot becomes available. 

The network 140 allows data and control information to 
be conveyed between the various infrastructure elements and 

10 operates in accordance with well-known network protocols. 
The fixed infrastructure may also include various other 
network entities 170 that may serve as the requesters of, or 
the destination for, location information for a given mobile 
unit. Examples of such infrastructure entities include, but are 

15 not limited to, consoles, telephone interconnect devices, and 
management terminals. 

The location processor 130 includes a digital computing 
device 131 and memory 132 used for the storage and execution 
of software programs. In practice, a commercially available 

20 computer workstation or a customized computer platform may 
be used to implement the location processor 130. In the 
preferred embodiment, the methods described hereinafter are 
implemented as software programs resident in the location 
processor 130 and/or mobile units 160. 

25 FIG. 2 is a block diagram of a location processor 200 

comprising a variance determiner -204, a differential variance 
determiner 206 coupled to the output of the variance 
determiner 204, and a location determiner 208 coupled to the 
output of the differential variance determiner 206. In the 

30 preferred embodiment, each block shown in FIG. 2 is 

implemented as one or more software-based algorithms stored 
in memory and executed by a digital computing device resident 
in the location processor 200. It is nonetheless understood 
that other, non-software based implementations may also be 

3 5 used. 
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The variance determiner 204 takes as input channel 
quality metrics. As described below, each channel quality 
metric uniquely corresponds to a signal received by a 
particular mobile unit or, in another embodiment, to a signal 
5 received by one of the fixed transceivers. The output of the 
variance determiner 204 is time of arrival variances having a 
one-to-one correspondence with the channel quality metrics. 
The time of arrival variances are indicative of the uncertainty 
as to exactly when each signal was received by the mobile unit 

10 (or fixed transceiver). 

The differential variance determiner 206 combines the 
time of arrival variances to produce time of arrival 
differential variances. As described below, the time of arrival 
differential variances are used to populate a time of arrival 

15 differential weighting matrix used in a WLS solution, as 
calculated by the location determiner 208. The location 
estimate output by the location determiner 208 represents an 
iteratively updated nominal location estimate judged to be 
sufficiently accurate when correction values used to update 

20 the nominal estimate have converged below minimum 
thresholds. 

In the preferred embodiment, the location processor of 
FIG. 2 cooperates with a mobile unit when determining a 
location for that mobile unit. FIG. 3 illustrates operation of a 

25 mobile unit when location information for the mobile unit is to 
be determined in conjunction with a location processor. At 
step 302, the mobile unit receives M signals from M different 
geographically-diverse transceivers (M £ 3). Each of the 
signals uniquely corresponds to one of the transceivers, and 

30 each signal comprises known data symbols transmitted by its 
corresponding transceiver. 

The delay between when the known data symbols were 
transmitted and when they were received serve as the basis 
for determining pseudo-ranges between the mobile unit and 

35 each of the transceivers. In practice, the delay calculated for 
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a given transceiver is actually an average performed over the 
delays measured for multiple receptions of the known data 
symbols. Referring to the example shown in FIG. 1, the mobile 
unit 160 is assumed to receive signals from three 
5 transceivers, labeled Xo, Xi and X2, each transceiver also 
having associated with it a pseudo-range Do, D1 and D2, 
respectively. 

As known in the art, differential pseudo-ranges (DPRs) 
can be formed and used as the basis for a WLS location 
10 solution. Continuing with the example of FIG. 1, DPRs can be 
defined as: 



15 



25 



D 02 =D Q -D 2 



Although the serving site is used as the reference in the above 
example, one of the neighbor sites can just as easily be used. 
In practice, the known data symbols are subject to errors, 
including time of transmission synchronization errors (TOTS 
20 errors). In order to account for such errors, the PR 

measurement between a given transceiver and the mobile unit 
is defined as: 



D i = D ± + e i 

and the DPRs are defined as: 



D 02 =D 0 -D 2 

30 

Where e± is the known TOTS error for the i'th transceiver. If 
PRs are calculated by the mobile, such TOTS errors for each 
transceiver can be conveyed to the mobile. 

Having received M signals from the transceivers, the 
35 mobile unit, at step 304, determines M channel quality metrics. 
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Preferably, an average over multiple receptions of known data 
symbols are used to derive each channel quality metric. In this 
manner, the effects of noisy signal quality measurements can 
be mitigated. The channel quality metrics take the form of 
5 carrier-to-interference-pius-noise ratios (C/l+N). Techniques 
for measuring such C/l+N ratios are taught in U.S. Patent No. 
5,440,582 entitled A METHOD AND APPARATUS FOR 
DETERMINING SIGNAL USABILITY, which patent is hereby 
incorporated by reference verbatim and with the same effect 

10 as though the same patent were fully and completely set forth 
herein. However, other indicia of signal quality (e.g., received 
signal strength) may be used in place of, or as supplement to, 
C/l+N. Regardless of how determined, each channel quality 
metric has a one-to-one correspondence to the known data 

15 symbols received from at least three sites (M > 3). In practice, 
one of the sites will be the current serving site for the mobile 
unit, although it is permissible that each of the at least three 
sites is a neighboring site to the current serving site. 

At step 306, the mobile optionally determines M RMS 

20 delay spreads, wherein each RMS delay spread uniquely 
corresponds to one of the M signals. As known in the art, 
received signals are often subject to distortion due to 
multipath interference and delays. An RMS delay spread 
expresses the standard deviation (or root-mean-square) value 

25 of the delay of reflections, weighted proportional to the energy 
in the reflected waves. In this manner, the RMS delay spread 
serves as an indicator of the effect of the current multipath 
propagation environment on known data symbols received by 
the mobile unit. Methods for determining RMS delay spreads 

30 are known in the art. 

At step 308, the channel quality metrics and the RMS 
delay spreads, if calculated, are sent by the mobile unit to the 
location processor. Assuming that multiple receptions of 
known data symbols are averaged together when calculating 
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each channel quality metric, the number of receptions of known 
data symbols used, N i$ is also transmitted. 

At step 310, responsive to the information sent in the 
previous step, the mobile unit optionally receives a location 
5 estimate for the mobile unit from the location processor. It is 
also possible that having executed steps 302-308, the mobile 
unit will not receive the location estimate. This would be the 
case, for example, when the location information is requested 
by some party other than the mobile unit, e.g., an 

10 infrastructure entity associated with the wireless 
communication system, or even another mobile unit. 

Referring now to FIG. 4, a flowchart is shown illustrating 
operation of a location processor when location information 
for a mobile unit is to be determined in conjunction with the 

15 location processor. At step 402, in response to a request to 
determine a location estimate for the given mobile unit, the 
location processor determines M channel quality metrics for 
the mobile unit. In the preferred embodiment, the channel 
quality metrics are "determined" by the location processor 

20 when they are received as inbound data from the mobile unit. 
In such a case, for each channel quality metric, the number of 
receptions of known data symbols used in the calculation of 
that channel quality metric is also received. However, as 
those having ordinary skill in the art will appreciate, it is 

25 possible that the location processor can determine suitable 
channel quality metrics itself based on other data received 
from the mobile unit. 

At step 404, the location processor may optionally 
receive information regarding at least one RMS delay spread 

30 measurement from the mobile unit. When provided, the RMS 
delay spread is preferably included along with the channel 
quality metrics and other information discussed in the 
previous step. 

At step 406, the location processor determines M time of 
35 arrival variances, one for each channel quality metric. To the 
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extent that known data symbols are used to derive the channel 
quality metrics, the time of arrival variances are indicative of 
the uncertainty in the measured propagation delay for each 
signal. In the preferred embodiment, the channel quality 
5 metrics are mapped into time of arrival variances according to 
predetermined relationships. FIG. 5 illustrates examples of 
such predetermined relationships. 

Each predetermined relationship illustrated in FIG. 5 is 
depicted by a curve 501-503. These predetermined 

10 relationships are derived using computer simulation and/or 
field testing to determine expected time measurement 
performance under a wide range of conditions. Using such 
curves 501-503, values of C/l+N (measured in decibels along 
the horizontal axis) are mapped into corresponding time of 

15 arrival variance values (measured in |isec 2 along the vertical 
axis). Value ranges for both the horizontal and vertical axes 
other than those shown are possible and are a matter of design 
choice. In practice, each predetermined relationship is 
implemented using stored tables or empirically-defined 

20 mathematical functions. Although a single family of curves 
501-503 is shown in FIG. 5, it is possible to derive many 
families of such curves. 

Assuming use of RMS delay spread measurements, each of 
the curves 501-503 is indexed according to a particular value 

25 of a received RMS delay spread. Each curve 501-503 

represents a range of RMS delay spreads centered around a 
single value. Thus, in the examples shown, a first 
predetermined relationship, represented by a first curve 501, 
is used when the received RMS delay spread is closest to 5.0 

30 \isec; a second predetermined relationship, represented by a 
second curve 502, is used when the received RMS delay spread 
is closest to 2.5 \isec; and a third predetermined relationship, 
represented by a third curve 503, is used when the received 
RMS delay spread is closest 0 y.sec. In the event that RMS 

35 delay spreads are not used, the single predetermined 
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relationship used would comprise either an average of the 
various curves or would be chosen as that curve corresponding 
to the most frequently measured RMS delay spread. 

Returning to FIG. 4, differential time of arrival variances 
5 are determined, at step 408, according to the equation: 

aj =fL +£1.^+1^ 
Thus, each differential time of arrival variance ^ is seen to 

°0i 

10 be the sum of a scaled time of arrival variance v 0 for a first 
site and a scaled time of arrival variance v ± for an Pth 
neighboring site to the first site (1 < i < M-1). The scaled time 
of arrival variance v Q for the first site is the time of arrival 
variance for the first site, 5.2, divided by the number of 

15 receptions of known data symbols, n q , sampled to determine 
the channel quality metric (and, subsequently, the time of 
arrival variance) for the first site. Likewise, the scaled time 
of arrival variance for the Pth site, v it is the time of arrival 
variance for the Pth site, ^2, divided by the number of 

20 receptions of known data symbols used, n ± . Although the first 
site is preferably the serving site for the mobile unit under 
consideration, it is equally possible for the first site to be one 
of the neighboring sites. 

Assuming that the mean DPR error is zero and that the 

25 known data symbol measurements at each site are 

uncorrelated, the elements of a time of arrival differential 
weighting matrix (or covariance matrix), v, can be determined 
according to the following equations: 

30 =u 0 +Vj for j = k; and 

v-k = v 0 for j * k. 

For a location determination involving a first site and M-1 
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neighbor sites, the time of arrival differential weighting 
matrix becomes: 









• ^0 










V = 














• • 



5 

recalling that d^-Vo+v/ Thus, the present invention 
provides a convenient method for estimating a covariance 
matrix for use in a WLS solution. 

Having determined a covariance matrix, the location 
10 processor, at step 410, solves for the WLS location solution 
given by the equation: 

YA D = VCAxy 

15 The procedures for determining the WLS location solution are 
further described with reference to FIG. 6 below. The result of 
determining the WLS location solution is a location estimate 
for the given mobile unit. Consequently, at step 412, the 
location estimate is optionally sent to a mobile unit or 

20 infrastructure entity that originally requested the location 
information. 

FIG. 6 illustrates a method for determining a location 
estimate using a WLS solution. At step 602, a pseudo-range 
differential correction vector A D and a differential direction 

25 cosine matrix c are determined based on a nominal location 
estimate, denoted by the coordinates (x n ,y n ). The initial 
nominal location estimate is selected to be within the 
boundaries of the first cell (preferably the serving cell), and is 
also selected based upon the type of cells being used. That is, 

30 if the first cell is "sectorized", as known in the art, the center 
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of mass of the serving sector is selected as the initial nominal 
location estimate. However, if the first cell is not sectorized, 
a good choice for an initial nominal location estimate is a 
point on a circle centered on, and equal to approximately half 
5 the radius of, the first cell. Furthermore, knowledge of which 
neighboring cells are to be used in the location determination 
may be used to further refine selection of such a point by 
possibiy selecting a point closest to the neighboring cells. 
Given the initial nominal location estimate and 
10 knowledge of both the transceiver locations and TOTS errors, 
e 0 and e if for the first site and i'th neighboring transceiver (1 
< i < M-1), respectively, nominal differential pseudo-ranges 
can be calculated according to the equation: 

15 n n ,oi = {£>n,o +e 0 )-(l>x, i i+e i ) 

wherein D n t0 is the known distance between the nominal 
location estimate and the serving site, and D nti is the known 
distance between the nominal location estimate and the 

20 transceiver for the i'th neighboring site. With the nominal 
DPRs, the correction values included in A D are calculated as 
the difference between the DPRs measured by the mobile (and 
conveyed to the location processor) and the nominal DPRs. 
Thus, for a location determination involving a first site and M- 

25 1 neighbor sites, A D is defined as: 
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The differential direction cosine matrix c includes 
differential direction cosines for each of the M-1 neighboring 
sites, and is defined as: 



c = 



c x ,n,01 Cyjifll 



c xji£>M-l c yji,0M-l_ 



Denoting the known location of the i'th neighboring site as 
( x i>Yi)> tne elements of c are defined as follows: 



10 



D n,0 



Cyjifli- 



_ y n - Y0 Yn ~ Yi 
Dn.0 D n.i 
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Having determined A D and c, a nominal location 
correction vector can be determined at step 604. The 
nominal location correction vector contains corrections for the 
nominal location estimate and is defined as: 



A xy = 



20 Using linear algebra to solve WLS location solution, the 

nominal location correction vector is given by the equation: 



A xy = \g T v V) 1 Q T v X A 



25 At step 606, having solved for the nominal location 

correction vector, the nominal location estimate is updated by 
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adding A^, i.e., + 5 x ,y n -h8 y ). With this new nominal 

location estimate, both the pseudo-range differential 
correction vector A D and the differential direction cosine 
matrix c are also recalculated to reflect the correction to the 
nominal location estimate. This process of calculating A D , c 
and A^; updating the nominal location estimate; and re- 
calculating A D , c and A^y can be iterated until a final solution 

is reached, as described below. 

At step 608, it is determined whether the nominal 
location correction vector is diverging. This is done by 
comparing the components of the current nominal location 
correction vector ($ x &2xmnt£yscumsnt) with the components of 
the previously calculated nominal location correction vector 
(S x , P zB V ±>us£y,pm V ±>iis)- ,n particular, divergence is deemed to 
have occurred when the following conditions are met: 

$x /ourtsnt > U ' ,pxev±>us 
Sy^uzxent > U ' §y,p2&v±>us 

u is an upper threshold having a value greater than unity and, 
in the preferred embodiment, has a value of 10. It is 
understood that, rather than a single upper threshold, separate 
thresholds, u x and u yf having different values may be used for 

the separate divergence conditions. If it is determined that 
divergence has occurred, further iterations of the WLS solution 
are discontinued, and no location estimate is determined. 

If the divergence conditions are not met, it is determined 
at step 612 whether a magnitude of the nominal location 
correction vector has converged below a lower threshold, as 
described by the equation: 
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The value of l is a matter of design choice, and in simulations 
has been set to a value of 8.05 meters (0.005 miles). If the 
nominal location correction vector has not converged below the 
lower threshold, processing returns to step 604 where the WLS 
5 solution is allowed to iterate another time. However, if the 
nominal location correction vector has converged, the location 
estimate is set equivalent to the current nominal location 
estimate at step 614. 

In light of the above, the present invention is seen to 

10 provide an improved method for estimating a covariance 

matrix, which in turn enables the computationally-efficient 
WLS solution to be used in a location determination context. In 
addition to the embodiment in which the location processor 
cooperates with a mobile unit, the present invention may also 

15 be beneficially employed in those cases where a location 

determination is desired using only the fixed infrastructure 
(transceivers and location processor) or only a mobile unit. 
These other embodiments are described with reference to FIGS. 
7-9. 

20 FIG. 7 illustrates an infrastructure-only approach to 

location determination in accordance with the present 
invention. In particular, the method illustrated in FIG. 7 would 
be suitably performed by the infrastructure illustrated in FIGS. 
1 and 2. At step 702, M different transceivers (M > 3) each 

25 receive a representation of a signal transmitted by a mobile 
unit, resulting in M different representations of the signal. 
That is, each transceiver receives the same signal, but subject 
to different error conditions (i.e., fading, interference, etc.) 
than the other transceivers. For example, referring to FIG. 1, 

30 transceivers Xo, ^ and X 2 receive a signal transmitted by the 
mobile unit 160. In the preferred embodiment, the signal 
transmitted by the mobile unit is a sequence of known data 
symbols essentially equivalent to those transmitted by the 
transceivers, as described above. 
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At step 704, each transceiver determines a channel 
quality metric for its corresponding representation in 
essentially the same manner as described above (FIG. 3, step 
304). Likewise, at step 706, an RMS delay spread may 
5 optionally be determined by each transceiver as described 
above (FIG. 3, step 306). Also, PRs are determined by each 
transceiver. At step 708, the channel quality metrics and, if 
avaiiabie, the RMS delay spreads, are sent by the transceivers 
to the location processor. Additionally, each transceiver 

10 conveys its PR, and the number of receptions of known data 
symbols used to determine the channel quality metrics to the 
location processor. 

Having received the information sent at step 708, the 
location processor executes steps 710-716 to determine a 

15 covariance matrix and WLS location solution. Because steps 

710-716 are essentially equivalent to steps 406-412 of FIG. 4, 
further description is not necessary. 

Another embodiment of the present invention, a mobile- 
only approach, is illustrated in FIGS. 8-9. FIG. 8 is a block 

20 diagram of a mobile unit suitable for a mobile-only approach to 
location determination. In particular, the mobile 800 
comprises a receiver 802, a variance determiner 804 coupled 
to the output of the receiver 802, a differential variance 
determiner 806 coupled to the output of the variance 

25 determiner 804, and a location determiner 808 coupled to the 
output of the differential variance determiner 806. The 
variance determiner 804, differential variance determiner 
806, and the location determiner 808 are structurally and 
functionally equivalent to the variance determiner 204, 

30 differential variance determiner 206, and the location 

determiner 208, respectively, described above with reference 
to FIG. 2. Also, the variance determiner 804, differential 
variance determiner 806, and the location determiner 808 are 
preferably implemented as one or more software-based 

35 algorithms stored in memory and executed by a digital 
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computing device resident in the mobile unit. For the purposes 
of the present invention, the sole difference between the 
location processor 200, and the mobile 800 is the inclusion of 
the receiver 802 in the mobile 800. At a minimum, the 
5 receiver 802 must be capable of receiving signals from the 
infrastructure transceivers and, in the preferred embodiment, 
must be capable of receiving known data symbols from 
different transceivers. Additionally, the receiver 802 includes 
the functionality, described above, required to extract channel 

10 quality metrics from the received signals. 

Operation of the mobile unit 800 is further illustrated in 
FIG. 9. In particular, steps 902-906 are equivalent to steps 
302-306, respectively, of FIG. 3. That is, the mobile unit 
receives signals from the transceivers and derives the 

15 information required for determination of a WLS location 

solution. Then, rather than transmitting such information to 
an infrastructure-based location processor, the mobile unit 
executes steps 908-914 which are equivalent to steps 406- 
412, respectively, of FIG. 4. 

20 As demonstrated above, the present invention 

beneficially provides a method and apparatus for estimating a 
covariance matrix for use in WLS location determinations. 
Using channel quality metrics as a basis, the present invention 
allows a covariance matrix to be efficiently estimated and 

25 employed in a WLS location solution. Because the WLS location 
solution is computationally efficient in comparison with 
techniques utilizing Kalman filters, and does not require the 
use of costly GPS receivers, the present invention may be 
implemented in either an infrastructure-based location 

30 processor or a mobile unit, or both in conjunction. Using 
computer simulations, it has been demonstrated that the 
present invention can determine the location for a mobile unit 
to within 125 meters (0.08 miles or 410.1 feet) of its actual 
location at least 67% of the time. 
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Although the present invention has been described in 
terms of the presently preferred embodiments, it will be 
readily apparent to those skilled in the art that various 
modifications can be made therein without departing from the 
5 spirit of the present invention or from the scope of the 
appended claims. 

We claim: 
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Claims 

1. In a wireless communication system comprising a fixed 
infrastructure in wireless communication with a plurality of 
mobile units, the fixed infrastructure comprising a location 
5 processor, a method for the location processor to estimate a 
location of a mobile unit of the plurality of mobile units, the 
method comprising steps of: 

determining M channel quality metrics corresponding to 
10 the mobile unit, wherein M > 3; 

determining M time of arrival variances based on the 
channel quality metrics; 

15 determining M-1 time of arrival differential variances 

based on the time of arrival variances; and 

determining a weighted least squares solution using the 
time of arrival differential variances to estimate the location 
20 of the mobile unit. 
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2. The method of claim 1, the step of determining the time 
of arrival variances further comprising a sub-step of: 

mapping carrier-to-interference-plus-noise ratios into 
5 unsealed time of arrival variances according to one of at least 
one predetermined relationship, wherein the one of the at least 
one predetermined relationship is indexed according to at least 
one RMS delay spread. 

10 3. The method of claim 1, further comprising the step of 
determining the time of arrival differential variances 
according to the equation: 

^ ^ a 2 - 

0i N 0 N ± U 1 

15 

wherein: o^. is an i'th time of arrival differential variance; 
v 0 is a first scaled time of arrival variance 
corresponding to a first site; 
v i is an i'th scaled time of arrival variance 
20 corresponding to an i'th of M-1 neighboring sites; 

d 2 Q is a first time of arrival variance corresponding 

to the first site; 

6^ is an i'th time of arrival variance corresponding 

to the i'th neighboring site; 

25 n 0 \s a number of receptions of known data 

symbols used in determining a first time of arrival 
measurement corresponding to the first site; 
Ni is a number of receptions of known data symbols 
used in determining an i'th time of arrival 

30 measurement corresponding to the i'th neighboring 

site; and 

the index i is in the range from 1 to M-1, inclusive. 
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4. The method of claim 3, further comprising the step of 
determining the time of arrival differential variances, wherein 
the first site is a serving site for the mobile unit. 

5 

5. The method of claim 3, the step of determining the 
weighted least squares solution further comprising sub-steps 
of: 

10 determining a pseudo-range differential correction 

vector A D and a differential direction cosine matrix c based 
on a nominal location estimate; 

calculating a nominal location correction vector A^y 
15 according to the equation: 



20 



wherein a time of arrival differential weighting matrix v is 
defined as: 

v 0 



V - 



01 
^0 

^0 



_v 0 v 0 v 0 ... a 2 ^ 
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updating the nominal location estimate, the pseudo-range 
differential correction vector, and the differential direction 
cosine matrix based on the nominal location correction vector; 
and 



iterating the previous steps of calculating and updating 
so long as the nominal location correction vector is not 
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diverging and so long as the nominal location correction vector 
has not converged beyond at least one lower threshold. 

6. The method of claim 5, the step of determining the 

5 weighted least squares solution further comprising a sub-step 
of: 

setting the location of the mobile unit equivalent to the 
nominal location estimate when the nominal location 
10 correction vector has converged beyond the at least one lower 
threshold. 

7. The method of claim 5, the step of determining the 
weighted least squares solution further comprising a sub-step 

15 of: 

discontinuing the iteration of the steps of calculating 
and updating when the nominal location correction vector has 
diverged beyond at least one upper threshold. 
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8. A location processor comprising: 

a digital computing device; 

5 a memory coupled to the digital computing device, the 

memory having stored thereon executable instructions that, 
when executed by the digital computing device, cause the 
digital computing device to comprise: 

a variance determiner, responsive to M channel 
1 o quality metrics corresponding to a mobile unit, that 

produces M time of arrival variances, wherein M > 3; 

a differential variance determiner, fesponsive to 
the time of arrival variances, that produces M-1 time of 
15 arrival differential variances; and 

a location determiner, responsive to the time of 
arrival differential variances, that calculates a weighted 
least squares solution to estimate a location of the 
20 mobile unit. 
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9. A mobile unit comprising: 



a receiver, responsive to M signals from M transceivers 
of a plurality of geographically-diverse fixed transceivers, 
5 that produces M channel quality metrics, wherein each of the 
signals and each of the channel quality metrics uniquely 
corresponds to one of the transceivers, and wherein M > 3; 



a location processor, coupled to the receiver, comprising: 
10 a digital computing device; 

a memory coupled to the digital computing device, 
the memory having stored thereon executable 
instructions that, when executed by the digital 
15 computing device, cause the digital computing device to 

comprise: 

a variance determiner, responsive to the channel 
quality metrics, that produces M time of arrival 
variances; 



20 



a differential variance determiner, responsive to 
the time of arrival variances, that produces M-1 time 
of arrival differential variances; and 



25 a location determiner, responsive to the time of 

arrival differential variances, that calculates a 
weighted least squares solution to estimate a location 
of the mobile unit. 
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10. A communication system comprising: 



at least M geographically diverse transceivers that each 
receive a corresponding one of M representations of a signal 
5 sent by a mobile unit, and that each determine a channel 
quality metric for the corresponding one of the 
representations, wherein M > 3; 

a location processor, coupled to the transceivers, 
10 comprising: 

a digital computing device; 

a memory coupled to the digital computing device, 
the memory having stored thereon executable 
instructions that, when executed by the digital 
computing device, cause the digital computing device to 
comprise: 

a variance determiner, responsive to the channel 
quality metrics, that produces M time of arrival 
variances; 



15 



20 



a differential variance determiner, responsive to 
the time of arrival variances, that produces M-1 time 
of arrival differential variances; and 

a location determiner, responsive to the time of 
arrival differential variances, that calculates a 
weighted least squares solution to estimate a location 
of the mobile unit. 
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